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_ , Abstract 

qj t We discuss how the Cholesky decomposition may be used to ascertain whether a critical point 

fJL^ ■ of the field theory scalar potential provides a stable vacuum configuration. We then use this 

| method to derive the stability conditions in a specific example. 

1 Introduction 

^p" 1 ' Understanding the structure of the vacuum is of central importance in quantum field theory. 

f^i |. Local minima of the effective quantum potential correspond to stable vacuum configurations of 

the system and physical states in the quantum theory are identified with fluctuations around a 
given minimum. 

In this note we shall introduce the Cholesky decomposition and present a method for using 
the Cholesky decomposition of the Hessian of the potential to ascertain the stability of a given 
vacuum configuration. Whilst approaches of this manner are familiar in some mathematical 
disciplines, it does not seem to be widely known in theoretical physics. 
00 ' This technique should be of value in the study of vacuum stability in models of interest in 

fN| , particle physics. In particular, the study of the vacuum structure of gauge mediated supersym- 

metry breaking models, where stability analysis is often complicated due to the number of scalar 
fields entering in the potential pQ. 

Firstly, we shall present a summary of the relevant notions from linear algebra and discuss 
their use in demonstrating vacuum stability. To illustrate this approach, in the latter half of the 
paper, we shall apply these techniques to a particular example for which vacuum stability has 
not yet been shown. Additional areas for potential application of this approach are suggested in 
the concluding remarks. 
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2 The Cholesky decomposition 

Given a general scalar potential V = V((j>x, • • ■ , (f> n ), & stable (degenerate) vacuum configuration 
of the system is a (semi)stable critical point of V. That is a set of parameters for which 

dV 

— =0 for i = l,---,n. (1) 

d(j>i 

The Hessian H of a potential V(<f>i, ■ ■ ■ , 4> n ) is defined to be the matrix of second derivatives of 
the potential with respect to each of the fields it depends upon: 

dv , . . dv 

ii I ; •. ; I- (2) 

dV dV 
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The Hessian is a useful object since the stability of a given critical point can be inferred 
from the properties of the matrix. Suppose that the point p is a critical point of V, then it is 
(semi)stable if and only if the Hessian of V at p is positive- (semi)definite [2]. If the Hessian H 
is singular at p, then using row operations H may be separated into a nonsingular symmetric 
submatrix M and a set of zero vectors. If the submatrix M is positive-definite, then H is 
positive-semidefinite. 

For a Hermitian or real symmetric matrix A the following three statements are equivalent 

®: 

(i) . A is positive-definite, i.e. 

(x,Ax)>0 (3) 

for all nonzero vectors x. 

(ii) . All of the eigenvalues of A are positive. 

(iii) . A can be uniquely decomposed as follows 

A = C^C (4) 
where C is an upper triangular matrix with real strictly positive diagonal entries. 



The factorisation of A into C and its conjugate in (iii) is referred to as the Cholesky decom- 
position of A and may be considered as the matrix analogue of the square root operation for 
scalar numbers. 

Additionally, the three statements below are also equivalent to each other [3]: 

(a) . A is positive-semidefinite, i.e. 

(x,Ax)>0 (5) 

for all nonzero vectors x. 

(b) . All of the eigenvalues of A are non-negative. 

(c) . There exists a (possibly singular) matrix B such that 

A = B^B. (6) 



Whilst one can ascertain whether a matrix is positive- (semi) definite by examining its eigenval- 
ues, finding analytically the eigenvalues of H evaluated at p becomes increasingly more difficult 
as the dimension of the parameter space grows. Demonstrating the existence of the Cholesky 
decomposition provides an alternative method to finding the eigenvalues of the Hessian and can 
often be more practical to compute. 

Since the Hessian is symmetric, if it is nonsingular at a given critical point p then it can 
be factorised as H = U'U, where U is an upper triangular matrix with symbolic entries. The 
requirement that the diagonal entries of U be real and positive provides the stability conditions 
for the vacuum configuration. Note that when the diagonal entries are real and positive, the 
matrix U is the Cholesky decomposition of H. If the Hessian is singular at p, then one must 
inspect the nonsingular submatrix M. Similarly, one may factorise M = Wll and the vacuum 
is semistable if the diagonal entries of U are real and positive. 

Note that the dimension of the singular subspace gives the number of massless states, some of 
which may be Nambu-Goldstone bosons, but others correspond to moduli ('flat-direction' fields 
in the supersymmetric case). 

For a given Hermitian or real symmetric matrix A with symbolic entries, there are algo- 
rithms [4], based on iterative elementary row operations, through which one may obtain the 
factorisation A = U^U. In the case where A is positive-definite this provides the Cholesky 
decomposition. Moreover, most standard mathematical computing packages can perform this 
operation via preprogrammed functions [5] . 
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3 A particular example 



We shall illustrate here the method outlined in Section [2] by deriving the stability conditions for 
the vacuum configuration of a model proposed by Chan and Tsou [5J, related to the Standard 
Model. The problem requires the analysis of a 10x10 Hessian and thus provides a sufficiently 
complicated example to demonstrate the utility of our method. 

The scalar content of the theory under consideration is as follows: one complex scalar field 
if transforming in the fundamental representation of su(2) and three complex scalar fields <j) 1 
(i = 1, 2, 3) transforming in the fundamental representation of su(3). The scalar potential of the 
model has the following form 

3 3 

v = - f i W \f\ 2 + \ w \f\ 4 - f i S j2J2\& 

a—l i—1 

+ a *(EEki 2 ) + ks E l^-cf 

\a— 1 i—1 / a, 6=1 

3 3 3 

a— 1 i—1 a—l 

The vector a is a global vector and is contracted against the multiplicity index i of the scalar 
fields (jf . A global symmetry of the model allows one to fix the orientation of a and henceforth 
we shall work with a = (1, 0, 0) T . Following [6J, we fix the symmetries of the theory and 
subsequently parametrise the scalar fields thus 

fr = 




X sin Si sin 7 e IX3 \ (8) 
Y cos 5 2 


X sin Si cos 7 e lX2 
Fsin 82 e %xi 
Z 

where the above variables are real and obey the following relationship 

X 2 +Y 2 + Z 2 = ( s 2 . (9) 
The quantities ( w and ( s are identified with the absolute lengths of the sets of scalar fields: 



c = v\W+\W+W\ 



(10) 
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Expanding the potential V in terms of equations {5} we obtain 



V — — fJ-wCw + ^wCw 

- [is (x 2 + Y 2 + Z 2 ^ + \ s (x 2 + Y 2 + Z 2 ' 

+ k s \x 4 + Y 4 + Z 4 
+ 2X 2 Y 2 sin 2 Sx cos 2 8 2 sin 2 7 
+ 2X 2 Y 2 sin 2 Sx sin 2 S 2 cos 2 7 
+ X 2 Y 2 sin 2 Si sin 27 sin 2<5 2 cos(x 2 — %i - X3) 
+ 2X 2 Z 2 sin 2 (5i cos 2 7 + 2F 2 Z 2 sin 2 <5 2 

+ VI C 2 (X 2 + F 2 + Z 2 ) - ^ 2 C 2 X 2 COS 2 (5! 



(11) 



For convenience the following quantities are defined 



A 1= X 2 - Y 2 , 



A 2 = Y 2 - Z 2 , 



(12) 



R = 



2 ' 



We consider the following critical point of the potential V, which we express in terms of Cs, Ai 
and A 2 : 

Xi = X2 = X3 = <*i = 5 2 = A 2 = 0, 
Ai =R, 

>2 _ 3 M ™ + (1 + 2R)u2g - 3^iC s 2 



C- 



6A W 

^2 _ 3/j. s - 3t/iC 2 + 
2k s + 6X S 

From equations (fT2")l and f| 13[) we obtain the relationships 



(13) 



1 + 2R 



Y = Z = Q 



1 - R 



(14) 



Since the variables X, Y and Z are real it follows that 

1 



(15) 



--<R<1. 

The Hessian H of the system is constructed as a symmetric matrix in the basis: 

{X, Y, Z, Cw, Si, S 2 , 7, tx, t 2 , r 3 }. 

The form of the Hessian at the critical point defined by the set of equations (TTBl is found to be 
the following 



H = 8- 



/ N 

^(2Y 2 n s + V2 C) 
K?- 



V" 











\ 



(16) 
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where N is a 4x4 matrix given by 



N 



/ X 2 (k s + X s ) XYX s 

XYX S Y 2 (k, + A.) 
XYX S Y 2 X S 



XFA S 

r 2 A s 
r 2 ( Ks + a s 

2 



2 



(17) 



The Hessian is positive-semidefinite if the submatrix TV is positive-semidefinite and the other 
entries of H are non-negative: 

< (2Y 2 k, + nC) , 

(18) 

< — re- 



using equations (|T4|) it can be seen that these entries are non-negative given the following 
conditions 

R - J 2 (19) 

K s > 0. 

Note that for certain values of the variables and parameters, the submatrix N is singular; for 
instance: 

R=l, -~, Cs = 0, or Cw = 0. 

For brevity we shall consider only the case where N is nonsingular, as the special cases where 
N is singular can be analysed using the same method and are computationally easier. 

Finding the eigenvalues of N analytically is difficult and has not been previously completed. 
Instead we shall utilise the method detailed in Section [2] to analyse the vacuum stability. The 
submatrix N may be factorised as N = U^U, where the form of U is given by 









( s X s y/(l-R) 



(,X s y/ (l-R) 



S +A s ) 



C 3 \/(l-fl)K 3 (K a +2A s ) Cs«sA a ^(l-fl) 



v' 3 ( K =+ A =) 






- v /3k s (k..-.\..J(.,-.,+2A 3 ~ 
C 3V / (l--R)K 3 (" a +3A s ) 
1 /3(k 3 +2A 3 ) 





-"2) 



CmQs^l + Aai^) 

2 v /k 3 (k 3 +A 3 )(k 3 + 2A 3 ) 

2^/k s (k s +2A s )(k s +3A 3 ) 
Cm i/4A„k| +(l2A s A ro -3i^i ! +2i/2i'i-i'|)K s -2A e ,i/; 
2yM^7+3A7) 



Recall that the quantities an d Cs are identified with lengths and thus arc defined to be 
non-negative. Whence we conclude that the diagonal entries of U are real and positive given the 
following further conditions 



Xs>- Y , 



X w > 



2v\X s + (3v 2 - 2v\V2 + v\)k s 
Ak s (k s + 3A S ) 



(20) 



In which case, by the argument detailed in Section^ the Hessian is positive-semidefinite. There- 
fore we have shown that the set of semistable vacua is nonempty and furthermore we have derived 
the provisions, (|T9|) and (|20|) . under which the vacuum configuration is stable. 



4 Concluding remarks 

We have outlined how the Cholesky decomposition may be applied to the Hessian of a system in 
order to ascertain the stability of a given critical point. In the example considered above it was 
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not practical to find the eigenvalues of the Hessian analytically and the Cholesky decomposition 
was key to obtaining the stability conditions for the model under consideration. 

Although techniques based upon the Cholesky decomposition are established in some math- 
ematical fields, they have not been widely used in theoretical physics for demonstrating vacuum 
stability or otherwise. One interesting application of the Cholesky decomposition was presented 
in [7], in the analysis of fermion masses in a wide class of effective low energy models emerging 
from intersecting D-brane configurations. 

Whilst the application of our technique is restricted in this note to the analysis of the Chan- 
Tsou model, it is evident from our calculation that the approach is quite general. This method 
is well suited for the study of models of interest in theoretical physics, many of which present 
a large number of scalar degrees of freedom, and it is expected that the procedure we have 
presented may find several further applications in particle theory. 

Supersymmetric theories are an immediate candidate for applying this technique as they 
generally involve many complex scalar fields, leading to complicated effective potentials. It 
is expected that the approach advocated above may be of use in stability analysis and the 
classification of flat directions in supersymmetric models [1] , [8] , [9] . 

A further opportunity for application is the stability analysis of various grand unified theories 
which feature, for phenomenological reasons, scalar fields in high dimensional representations of 
the gauge group |10] . Whilst work has been undertaken to identify stable minima in such models 
[llj . there remain open problems which might be approached using the method presented here. 
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